Manifestation of a gap due to the exchange energy in a spinor condensate 
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We investigate the dynamic response of population transfer between two components of a finite 
temperature spinor Bose condensed gas to a time-dependent coupling potential. Comparison be- 
tween results obtained in the Bogoliubov- Popov approximation (BPA) and in the generalized random 
phase approximation (GRPA) shows noticeable discrepancies. In particular, the inter-component 
current response function calculated in the GRPA displays a gapped spectrum due to the exchange 
interaction energy whereas the corresponding density response function is gapless. We argue that 
the GRPA is better since, contrary to the BPA, it preserves the SU(2) symmetry and the /-sum rule 
associated to the spinor gas. In order to validate the approximation, we propose an experimental 
setup that allows the observation of the predicted gap. 
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I. INTRODUCTION 

A great achievement in the discovery of condensation 
of Bose gas has been the experimental observation of the 
gapless and phonon-like nature of the collective excita- 
tions [1-3]. The experimental measurement of the en- 
ergy spectrum as a function of momentum has confirmed 
the validity of the Bogoliubov theory However, well 
understood at zero temperature, the finite temperature 
extension of this theory is still controversial. The main 
difficulty is to find an approximation both conserving and 
gapless @. The generalized random phase approxima- 
tion (GRPA) is such an approach that has been intro- 
duced first in the context of a quantum plasma in order 
to explain the collective phenomena [6j. Its application 
to Bose condensation has been successful to determine 
the sound velocity, the low-lying excitation spectrum ob- 
served in a trapped Bose gas [7| and the critical velocity 
@. Nevertheless, this theory suffers from the apparent 
contradiction that its building blocks contain propaga- 
tors describing the particle motion in terms of the gapped 
and parabolic Hartree-Fock (HF) mean field dispersion 
relation 0, H| . This paradox is resolved if the theory is 
reinterpreted adequately by making the distinction be- 
tween collective excitations that result from oscillation of 
the condensate and quasi-particles or thermal excitations 
that compose the normal part of the gas @. Generally, 
in other fluids the two kinds of excitation exist (e.g. the 
plasmons and the electrons in a plasma). Therefore, no 
reason justifies that they should be identical in a Bose 
condensed gas as predicted from the Bogoliubov theory. 

The only difference is that, in other fluids, the pole 
of the one particle Green function corresponds to the 
single particle spectrum while the pole of the suscepti- 
bility function corresponds to the collective excitations. 
In a Bose gas, below the condensation point, the two 
poles mix to form new poles that are identical for the 
two functions [5, 7]. In the non-conserving Bogoliubov 
theory, only one branch of excitation appears as a pole 
which corresponds to the quasi-particle resulting from 
hybridization of quantum states. In the GRPA, however, 



the two branches of the susceptibility function and the 
branch of the one particle Green function mix to form 
three common branches of excitations, one of which cor- 
responds to the phonon-like excitation Q- In that case, 
it is not possible to associate a single quasi-particle to the 
one particle Green function with three branches of excita- 
tion as poles. For these reasons, the idea that the normal 
gas is composed of Bogoliubov-like quasi-particles seems 
too simplistic and must be reconsidered in the frame- 
work of a conserving approach. Basically, the GRPA 
tells that the excitation spectrum of a condensed parti- 
cle is phonon-like while the one for the thermal particle 
is parabolic- like with a gap d, |]| . The whole difficulty 
is how to perturb the Bose gas in order to selectively 
address these excitations. 

In this paper, we aim at verifying the validity of these 
hypothetical thermal excitations by proposing an exper- 
iment on a spinor Bose gas revealing the existence of the 
gap. After reminding about the GRPA theory in section 
2, we study in section 3 the inter-component transition 
between two internal levels of the Bose gas subjected to 
an external coupling field [lfjl - fl3l |. In particular, we are 
interested in determining the energy required to transfer 
a thermal atom from a Bose condensed gas towards an- 
other internal sub-level. According to the GRPA, such an 
atom becomes distinguishable from the condensed atoms 
and consequently releases a fixed amount of energy due 
to the Fock exchange interaction. An analog experi- 
ment is encountered in a superconductor-metal junction 
in which the measure of the electrical current as function 
of the voltage allows the determination of the BCS gap 
Both experimental [ll| and theoretical studies |1(J | 
have been carried out on transitions between two inter- 
nal levels in order to probe the inter-component density 
response function. But an external coupling has the ef- 
fect to rotate the whole gas in the spinor space at no 
energy cost in the long wavelength limit and thus this 
function is gapless. For this reason, we investigate in- 
stead the inter-component current response function in 
order to be able to address separately the thermal atoms 
from the condensed ones and to display the gap. For com- 



2 



parison, a similar study is carried out in the gapless but 
non conserving Bogoliubov-Popov approximation (BPA) 
[!, [j| . This approach extends the Bogoliubov theory to 
strongly depleted gas and leads to different results. In 
section 4, we generalize the approach to the case of real 
experiments with a trap while section 5 ends with the 
conclusion. 



II. THE GRPA FORMALISM 

To start with, we consider a spinor Bose condensed gas 
of atoms with mass m in a volume f2 populated only in 
two sub-levels chosen among the hyperfine structure of 
the atom and labeled by a — 1, 2. The scattering lengths 
between atoms of any sub-level are assumed to be equal 
to a (in 87 Kb their relative difference is about 10 -2 |12| ) 
and define the coupling constant g — Aita/m {h = 1). 
Defining the annihilation operators c Qi k where k is the 
momentum and £k = k 2 /2m, the Hamiltonian is written 
as: 

H = ^2 £kC l,k C Q,k + ^ ^2 C l,k C lj,k' C Q,k+q C ^,k'-q (!) 



a.k 



CK,/3,k,k',q 



The eventual energy shift between the two sub-levels 
is removed by choosing an appropriate representation 
of c Q ,k- We shall study the dynamic evolution of the 

excitation operator p^qW = e%Ht <^ a k c ^,k+ q e ~ lHt in 
terms of which we express the space Fourier transform 
of the densities and the currents, respectively pq'^(t) = 

EkP£,'J(*) and J«."(t) = Ek<*S(*)(k + q/2)/m. These 
are associated to the U(l) number conservation and 
SU(2) spinor symmetries of H. We first consider the case 
of a uniform condensed gas initially stable populated in 
the sub-level 1 such that (Pk q) = ^a,i^a,/9^q,o^i,k- (• • • ) 
denotes the statistical ensemble average Tr(<r . . . ) where 
a is the density matrix, n^k the momentum density in 
the mode k and n c = n^o the condensate density. 

Let us introduce a coupling between the two spinor 
components through the time dependent Hamiltonian 
H v {t) = n P 1 ^V^l,e- i ^ +m)t + c.c. where V™ u is a con- 
stant and iO ensures the adiabatic condition. The change 
in the frequency mode u> of the inter-component density 
dPcibj from its equilibrium value (p q 2 ) = induced by 
this potential is related to the inter-component suscepti- 
bility function through: 



Sp, 



(2) 



Formal linear response theory allows to find the formula 
0: 

f°° + 
Xi2(q,«) = ^ / dt e j( ^ 0)t (K 2t (0),p q 2 (i)]) (3) 



Although the initial state is not the 517(2) symmetric 
equilibrium state, the linear response formalism is valid 



in a much more general context of a stable state such 
that [H, a] = 0. Similarly to the intra-component case, 
provided that this function is causal and the r£ a ,k' s are 
isotropic, its imaginary part obeys the following f-sum 
rule [101] : 



f°° q 2 

/ dwwlmxi2(q, = -7T— 
J-oo 2m 



Tlx 



(4) 



This function is also connected to the longitudinal cur- 
rent response function %| 2 (q, w) through the exact rela- 
tion: 



Xi2(q,w) = 



rii 



uj + iO 2m(uj + i0y 



+ 



q 2 xl 2 (q^) 

(uj + iO) 2 



(•5) 



This last function is defined by replacing p q 2 (i) by 
q.J q 2 (i)/|q| in Q. The /-sum rule and Eq.© result 
from the conservation law d t p^(t) = zq.J q 2 (i) associated 
to the SU{2) symmetry. 

The GRPA approach is based on the assumption that 
the off-diagonal part of the excitation operator (q ^ 0) 
contributes less than the diagonal one (q = 0) because 
the operators involved do not oscillate in phase but rather 
randomly Q . As a consequence, in the Heisenberg equa- 
tion for p k 2 qW> we neglect all contributions quadratic 
in the off-diagonal operator in comparison with products 
involving off and diagonal ones. In this way, we obtain 
an equation of motion of a similar form to the intra- 
component case [1, [l3| : 



(id t 



f 2,k+q 



HF 
l,k 



)Pk 2 q W=5n likPq 2 W 



(6) 



where = e k + <mi(l + is the HF energy felt by 
the atom in sub-level a. 



III. RESULTS FOR THE HOMOGENEOUS GAS 

The operatorial integral equation ^ is solved exactly 
using methods developed in [9|. Plugging the solution 
into (|3|) , the susceptibility function is found to be [7j| : 



Xi2(q,w) = xo,i2(q,w)/(l -3Xo,i2(q,w)) 



where 



Xa 



ia(q,w) = 



ni,k 



iO 



C HF 
-l.k 



HF 
fc 2,k+q 



(7) 



(8) 



As shown in [lJJ], this function is gapless since the pole 
uj = corresponds to a rotation in the SU(2) space. 
Using the Kramers-Kronig relations [4|, we check that 
Eq.© fulfills the /-sum rule Also, by introducing a 
vector potential as an external perturbation, the longitu- 
dinal current response function can be calculated as well 
in this approximation and the result verifies the property 
©. Combination of ©, © and jSJ in the limit of long 
wavelength allows to obtain: 



xi 2 (q = o, w ) 



^ k 2e kni,k/(3m) 



iO 



■gm 



(9) 
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We find a non zero pole — gn\ that corresponds to a 
gap due to the Fock exchange interaction and that is 
interpreted as the transition potential energy for any 
transferred thermal atom. The fluctuations of the inter- 
component longitudinal current induce a relative motion 
between the normal fluid and the superfluid in the spin 
space. Note that, contrary to the intra-component case, 
the longitudinal current response function Eq.© does 
not have the same pole as the density response function 
Eq. (JTJ) and is identical to the transverse current response 
function when q > 0. Eq.© is the main result of this 
paper and we shall show how this gap can be observed 
in experiments with a trap. 

In contrast, the use of the BPA leads to drastically 
different results. It corresponds to the Hartree-Fock- 
Bogoliubov approximation to the Hamiltonian in which 
the term quadratic in the product of anomalous aver- 
age is also neglected The chemical potential is given 
by // = g{2n\ — n c ), the quantum operators describing 
the various components evolve according to c±.k(t) = 

e-^tVU^Skfi + u+,ke- ie **&i, k + +u- M e K ^b\^) and 
C2,k(i) = e _i ^ I+ek ^*C2.k- &i,k is the annihilation operator 
associated to the quasi-particle, e k = \j2gn c t\ L + e k and 
M ±l k = ±((ek + S"c)/2e k ± 1/2) 1 / 2 . As a result, using 
(|3]), we obtain: 



xf 2 p (q^) = 



Xl2 (°> w ) 



iO 



£ 

±.k 



±,k 
B 



±, k v i, k ±,_^ 



iO±e* 



fc k±q 



2e k u±, k K,k + *±,-) 



3 m 



(11) 



where nf k = (b\ k b 1]l ). Both approaches GRPA and 
BPA are gapless and their susceptibilities defined respec- 
tively by © and (fTO)) have a pole at u = e q for zero 
temperature in agreement with the spectrum obtained 
from the linearization of the coupled Gross-Pitaevskii 
equations [16J. Besides this pole, we note from Eq. flTOl 
[TT]) that transitions in BPA result from scattering of a 
phonon into a normal atom releasing energy e k — e k+q 
and from creations of both a phonon and a normal atom 
absorbing the energy e k + e k _ q . As a consequence, for 
q = 0, the transition frequency spectrum is continuous 
with a line- width estimated to gri\. In GRPA, however, 
only a process that releases the constant gap energy hap- 
pens without any line-width. 

Therefore, the two theories display noticeable differ- 
ences but strong theoretical arguments play in favor of 
GRPA. Due to the SU(2) breaking of symmetry, the BPA 
violates the conservation laws and thus does not obey the 
/-sum rule nor Eq.©. In particular, for q = and lu = 0, 
the action of the external coupling is to rotate the whole 
gas in the spin space while it rotates only the conden- 
sate in the BPA. Physically, this violation can be under- 
stood by noticing that a phonon has not a well defined 
atom number ([&i.kj pj^o] ^ ^i,k) an d therefore cannot 
be transformed into an atom alone during a transition. 



For these reasons, the GRPA is a better approximation 
that is justified if the predicted gap is observed in exper- 
iments. 



IV. RESULTS FOR THE TRAPPED GAS 

The considerations so far obtained are straightfor- 
wardly extended to a trapped gas. We start with a con- 
densed gas initially populated in the sub-level 1 confined 
in a parabolic potential V n (r) — J2i= x y z \ m ^ 2 r 2 . 1 R 
this new situation, the density profile n\(v) depends on 
the position and is determined in good approximation 
from the mean field HF equations in the Thomas-Fermi 
limit 0| : the condensate density n c (r) is determined from 
the Gross-Pitaevskii equation in absence of the kinetic 
term /i = y n (r) + g(2rii(r) — n c (r)); the thermal atoms 
density is determined from the semi-classical expression: 
ni)k (r) = l/[exp(j3(eg F (r) + V n (r) - fi)) - 1] where 
/3 = 1/ksT is the inverse temperature. Would we have 
applied the BPA, we should have rather used nf k (r) = 
l/(exp(/3e k (r)) — 1) in the condensate region. We assume 
an isotropic trap with frequency Ui = uj = 2tt x 100Hz, 
scattering length a\JrnU5 = 10~ 3 and fi — 2tt x 2kHz. 
In the second sub-level, we consider both cases of a 
trapped and a fiat potential V 22 (r) = V n (r),0. In or- 
der to induce a current in the trap and transfer atoms 
to the second sub-level, we apply the external coupling 
during a finite time interval and 



V 12 {v,t) - r z 

beginning at t = where F z refers to a force constant 
along the z axis in the real space. This coupling can 
be produced either by a rotating magnetic field gradient 
B = B'r z (cos(u>t)t),sm(ujot),0) that interacts with the 
gas with its spin orientation along the z-axis or by 
a two photon Raman scattering between two multiplets 
(e.g. F = 1 and F = 2 in S7 Rb) creating a stationary 
optical field with a node centered at the minimum of the 
trap [l3| . 

In presence of inhomogeneity, the density response has 
to be generalized into 

<5p 12 ( ri , t) = f dt'f d 3 r 2 X u(ri,v 2 ,t - t')V 12 (r 2l t')(12) 
Jo Jn 

The total number of output atoms in the second level 
A2 = J n (i 3 rn2(r, t) is determined through the transfer 
rate equation: 



dN 2 
dt 



= i d^^V 2 
'n 



(r l! t)5/ 2 (r 1 ,t) + c.c. (13) 



For t ^> l/wo, the time integral can be ex- 
tended to infinity in (fT2|) and Eq.(fT3|) becomes time- 
independent. Combining Eqs. (11211131) . rewriting ri = 
r + r'/2 and r 2 = r — r'/2 and using xi 2 (r, q, w) = 
/^(//r'e-'^'+^xijfri,^^), the transfer rate 
is expressed up to the second order in the coupling in 
terms of the inhomogeneous susceptibility Xi2(r, q, lo). 
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FIG. 1: In (a), schematic representation of the local poten- 
tial energy (full line) in a trap felt by thermal atoms for the 
trapped sub-level 1 (left) and 2 (right). The difference ac- 
counts for the transition energy. In (b), corresponding rate of 
output atoms in the second sub-level in GRPA (black) and in 
BPA (grey) for T/T c = 0.56. 



This last function is calculated in the local density ap- 
proximation by replacing the density and the frequency 
in the homogeneous expressions ([5]) and ([§]) by their local 
forms m(r) and w(r) = V n (r) - V 22 (r) - w [3- This 
procedure gives successively: 



1 dN 2 
~F 2 ~dT 



2tt d 6 r 



= -2nJ tfrlm 
'm(r) 



d 2 \ 

+ -f - I Imxi2(r,q,w )| q ^. 
r>i(r) , 1 



4m 



uj(r)+i0 (uj(r)+iQ) 2 
Ek e k»i,k(r)/(3TO) 
u)(r) + i0 + gni(r) 



(14) 



Eq.(fT4"| is expected to describe accurately the transfer 
rate provided lJ and t <C uJ -1 so that the effects 

of discretization in the excitation spectrum can be safely 
ignored. These effects also diminish due to the higher 
density states at finite temperature. Also, the linear re- 
sponse theory is valid provided that the typical time du- 
ration t is much lower than the inverse of the typical Rabi 
frequency estimated to F z Rtf where Rtf is the conden- 
sate size. This implies the upper bound F z <C w/ Rtf ~ 
w 2 -\/m/2/i. Thus, in GRPA, the total number of trans- 
ferred atom results from a local rotation of the gas in 
the spin space in regions where ljq = V u (r) — V 22 (r) 
and from an inter-component current of thermal atoms 
releasing the exchange Fock interaction energy in regions 
where w = gni(r) + V n (r) — V 22 (r). 

In the case of V 22 (r) = V^ 11 (r), rotation in the spin 
space occurs only for zero frequency so that a non zero 
frequency rate provides a direct evidence of the gap 
cj = gn\{r). This released energy depends mainly on the 
local condensate density (see FigQJ,). Thus the transfer 
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FIG. 2: In (a), same as Fig.l but for a transition to an un- 
trapped sublevel. Note the constant gap in the condensate 
region responsible for the peak. Corresponding rate of out- 
put atoms for T/T c = 0.56 in (b) and T/T c = 0.85 in (c). 



R = m(ZJ/F z ) 2 dN 2 /dt calculated numerically from f| 14[) 
and shown in FigJTjD as a function of ujq displays an in- 
homogeneous broadening. This rate decreases with fre- 
quency whereas the density of thermal atoms is higher in 
the edge of the condensate than in the center of the trap. 
For comparison, we plot the transition rate obtained in 
BPA by inserting Eq. fTUf in Eq. (fT4"]) where contrary to 
GRPA energy absorption is possible when uio < 0. The 
reduced transition rate results from the broad spectrum 
displayed in Eq.(|TD|). 

FigUl shows the same curves as FigQ]but for the case 
V 22 (r) — 0. The output atoms are not trapped in or- 
der to diminish the dispersion due to the inhomogeneity. 
As a result, a resonant transition occurs for a frequency 
close to loq — gni(r — 0). However the transition due 
to the spinor rotation becomes frequency dependent and 
leads to a background noise. The broadening of the res- 
onance is caused by the variation of the thermal density 
in the condensate region. The line-width is estimated to 
Aoj ~ g(mfcsT/27r) 3 / 2 and increases with the tempera- 
ture. Such a broader peak is observed in Figl^t but at the 
expense of a more intense background noise. This noise 
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is proportional to the local gas inertia r 2 z ni(r) about the 
z-axis (see Eq. fLll) ) and is more important for stronger 
depletion in the region of high frequency. It considerably 
reduces if we choose uj z 3> 5J without altering the validity 
of the local density approximation (u) z <wo). 

V. CONCLUSIONS 

In conclusion, we analyze theoretically the atom tran- 
sition rate between two sub-levels of a spinor condensate 



caused by a time dependent external coupling. Signif- 
icant differences are reported between the GRPA and 
BPA models. An experimental setup is suggested in or- 
der to validate at most one approach and can serve to 
understand the true nature of the quasi-particle disper- 
sion relation. 
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